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Abstract. Structure of mean pulsar radiation patterns is discussed within the nested-cones and patchy beam 
models. Observational predictions of both these models are analyzed and compared with available data on pulsar 
waveforms. It is argued that observational properties of pulsar waveforms are highly consistent with the nested- 
cone model and, in general, inconsistent with the patchy beam model. 
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1. Introduction 

One of the important questions in pulsar research is what 
the overall structure of the mean pulsar beam is and how 
this structure is related to highly fluctuating instanta- 
neous pulsar radiation. It is difficult to reveal this struc- 
ture as pulsar observations represent one-dimensional cuts 
through two-dimensional beams. However, some indirect 
methods have been applied in an attempt to resolve this 
problem and two major models of pulsar beams have 



emerged from this work. Rankin (1993), Gil et aL (1993) 
and Kramer et al. (1994) calculated the opening angles p 
of emission corresponding to a pulse width W measured at 
10 and 50 percent of the maximum intensity. As a result, 
they obtained a binomial distribution of these angles, that 
is, for a given period P one of the two preferred values was 
possible, following however a general P~ x l 2 dependence. 
Such distribution is most naturally interpreted as an indi- 
cation of two nested cones in the structure of mean pulsar 
beams. This interpretation is called a conal model of pul- 
sar beams. An alternati ve model postulates that th e mean 
pulsar beam is patchy ( Lyne fc Manchester, 1988 , LM88 
hreafter) , with different components randoml y distributed 
within an almost circular "windo w function" ( Manchester, 
19951 ; |Han fc Manchester, 200"l| ). Such a model is appar 



ently inconsistent with the binomial distribution of the 
opening angles inferred from measured pulse widths. In 
fact, unless putative patches are distributed along nested- 
circula 



paLLt 



lems, the drslubuliuu uf euiiespuudiiig upeu- 
ing angles shuuld be (fui any given peiiud) laiiduui lalhei 

than binomial. 

Recently, [Vlitra fc Deshpande (199S| , MD99 hereafter) 
attempted to test both these rival models. They dis- 
tributed locations of the profile components (measured 
as the peak-to-peak separation of the outer conal compo- 



nents in complex profile pulsars) on one quadrant of the 
beam represented on the common normalised scale (with 
x-axis and y-axis representing longitudes ip and impact an- 
gles (3, respectively; refer to Fig. 4 of MD99). They found 
that most of the peak intensity points are concentrated in 
narrow sections of the beam. This feature is a strong indi- 
cation of the conal structure of the pulsar beam. It can be 
argued that if, indeed, the pulsar beams are patchy, then 
in such a case there is a high probability for the beam to 
be uniformly filled with peak intensity locations. Further 
in their analysis, they excluded the so-called conal sin- 
gle and conal double profiles (Rankin, 1983), which are 
thought to be exclusively grazing cuts of the line-of-sight 
at the beam boundary. They found that such exclusion 
in their sample led to the absence of points at high im- 
pact angles (3, which is perfectly consistent with the nested 
cone model and inconsistent with the patchy beam model. 
In fact, within the patchy beam model there is no reason 
why the single and double profiles occur exclusively at 
high impact angles. Moreover, the midpoint of single and 
double profiles usually coincides with the fiducial phase, 
at which the multifrequency profiles align after being cor- 
rected for cold plasma dispersive delays. This property is 
natural within the conal model and inconsistent (in gen- 
eral) with the patchy beam model, since patches would 
have to be placed symmetrically with respect to the fidu- 
cial plane, containing both magnetic m and spin O axes 
(Fig. |l|). As a rgued by MD and independently by Gil fc 
Scndyk (200C , GSOO hereafter) the "mean" average pulsar 
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beam consists of up to three nested cones, centered on the 
global magnetic dipole axis. 

Observations of single pulses in strong pulsars show 
that longitudes of subpulses are weakly dependent on fre- 
quency as compar ed with longitudes of c orresponding pro - 
file components ( Izvekova et al., 1993| ; Gil et al., 2002 ). 
Within the conal model, the longitudes of profile com- 
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Fig. 1. The beam emission patterns for both the conal model (illustrated on the right-hand side), and for the patchy 
beam model (illustrated on the left-hand side) are presented in this composite picture. The geometry of observation 
is determined in both models by the inclination angle a between the magnetic m and the spin SI axes and the impact 
angle j3 of the closest approach of the line-of-sight (observer O) to the magnetic axis m. The fiducial plane containing 
the spin axis SI, the magnetic axis m and the observer O defines the fiducial longitude = 0° which is common for 
both models. All pulse phases (ip s - longitude of subpulse peak, <p p - longitude of profile component peak) are measured 
from this fiducial longitude, to the left from it in the patchy model and to the right from it in the conal model. In this 
paper we ignore the dispersive delays which are the only cause of frequency variation at or near the fiducial plane. 
The frequency-dependent position of any observed feature in the beam pattern is described by two angles: the opening 
angle p(y) between the m axis and the line-of-sight (1-of-s) and by the azimuthal angle a between the fiducial plane 
and the plane of dipolar field lines associated with a particular feature. Within the conal model the subpulse emission 
is associated with the frequency dependent subpulse spots S(v), which move circumpherentially to form a cone C(y). 
Thus, the profile components are associated with the frequency dependent cones in the conal model, while within the 
patchy model the subpulse-associated spots SP(v) occupy the patchy areas PA(v) related to the profile components. 
The longitudes of subpulse peaks tp s (v) and profile peaks f P {v) are marked, with v\ > V2, in both models. 
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ponents arc determined by the intersection of the line of 
sight trajectory with the frequency-dependent cones of the 
maximum average intensity, while the longitudes of sub- 
pulses are determined by the intersection of the line-of- 
sight trajectory with subpulse-associated emission beams, 
which move across the average cones as frequency changes 
(Gil & Krawczyk, 1996). We demonstrate in this paper 



that the different frequency dependence of subpulse and 
profile component longitudes is a natural property of the 
conal model, and that both subpulses and profile compo- 
nents should demonstrate the same frequency dependence 
of their longitudes within the patchy model. We present 
both general qualitative arguments and detailed quantita- 
tive model calculations to support the above statements. 
For better understanding of our argum ents, this paper 
should be studied along with the paper by |Gil et al. (2002 , 
GGGK hereafter), in which details of frequency depen- 
dence of emission patterns in PSR B0329+54 are discussed 
within the conal model of pulsar beams. 

2. Geometry of pulsar radio beams 

It is widely believed that narrow-band pulsar emission is 



relativ stically beamed tangentially to dipolar magnetic 



field lines. Thus, the emission beaming geometry can be 
described by an opening angle p = 1°.24 s r~J 2 P~ 1 l' 1 , 
where r§ = r{y)/R is the normalized emission altitude (in 
units of stellar radius R = 10 6 cm). The mapping parame- 
ter 0<s = e?/r p <lis determined by the locus of dipolar 
field lines on the polar cap (s = at the pole and s = 1 
at the polar cap edge), where d is the distance from the 
magnetic axis m to the field line on the polar cap corre- 
sponding to a certain detail of the pulse profile (peak of 
subpulse or profile component), and r p = 1.4-10 4 F -1 / 2 cm 
is the canonical polar cap radius. According to the gen- 
erally - accepted concept of the radius-to-frequency map- 
ping, higher frequencies are emitted at lower altit udes r(y) 
than lower frequencies. Kijak fc Gil (1997, 1998) found a 



semi-empirical formula describing the altitude of emission 
region corresponding to a given frequency i/qhz (in GHz) 
which reads r$ w 50 • Vq^ 1 ■ t$ ' ■ P - 33 , where r 6 is the 
pulsar characteristic age in million years and P is the pul- 
sar period. This formula for emission altitudes is used in 
our model calculations. 

To perform geometrical calculations of the radiation 
pattern one has to adopt a model of instantaneous en- 
ergy distribution on the polar cap. Any specific inten- 
sity distribution can be transferred from the polar cap 
along dipolar field lines to the emission region, and then 
along straight lines (following the opening angles p = 
l°.24(d/r p )r6 /2 F- 1 /2) to 

a given observer specified by the 
inclination and impact angles (a, (3). We assume that at 
any instant the polar cap is populated by a number of 
features with a characteristic dimension V, delivering to 
the magnetosphere corresponding plasma columns flow- 
ing along separate bundles of dipolar magnetic field lines. 
Each feature can be modelled by the Gaussian inten- 
sity distribution. Since the elementary pulsar radiation 



is relativistically beamed along the magnetic field lines, 
we can transform the feature-associated intensity pattern 
to the radio emission region and obtain the subpulse in- 
tensity I observed at the longitude ip in the form = 
cxp(— nlf(ip) /T> 2 ) where the subscript i lebels different fea- 
tures and li(tp) = d 2 (ip) + df — 2d((p)di ■ cos[a(ip) — (Tin] ■ 
Each features is located at the position (Fig. Q) described 
by the polar co-ordinates di (distance from the pole) 
and <Ji n = Oio + nD r (magnetic azimuth angle), where 
the subscript o refers to the initial position correspond- 
ing to the first pulse, n is the sequential pulse num- 
ber and D r is the drift rate. The instantaneous emis- 
sion of the rt-th pulse is described by — 52i=ili(<p)i 
where the sum includes a number of adjacent features 
contributing significantly to the observed intensity, de- 
pends strongly on the inclination a and the impact (3 
angles, determining the cut of the line-of-sight tracjec- 
tory across the beam (Fig. []]). In fact, the running po- 
lar co-ordinates along the line-of-sight trajectory can be 
expressed in the form d(ip) = [p(ip)/l° .24}r p r 6 1 / 2 p 1 / 2 
and a(<p) = atan { c ^+py-c^ac^pl v ) }> an d p(tp) = 
2asin {sin 2 (y/2) sin a sin (a + 0) + sin 2 (/3/2)} , ( jGil et al 



1984). The average pulse profile is therefore I{ip) 
jf'E n f =1 I n ((p), where N is the number of averaged single 
pulses. The geometry of pulsar radiation described above 
(for more details see |Gil fc Krawczyk, 1996| , |l~997| , GGGK) 
can be applied to both conal and patchy beam models 
(Fig-0). 

We adopt a model of a pulsar beam in which the in- 
stantenous subpulse emission corresponds to a number of 
isolated subpulse beams, while the average emission re- 
flects the conal structure resulting from circumferential 



motion of subpulse beams (Rudcrman fc Sutherland, 1975 ; 
|Gil fc Krawczyk, 1996| , |l997| ; |Gil et al., 2002|) . Therefore, 



the longitudes of subpulse peaks correspond to phases of 
interception of the subpulse beams by the line-of-sight tra- 
jectory, while the longitudes of profile component peaks 
are determined by the intersection of the line-of-sight tra- 
jectory with the average cones. The important point is 
that within such model the subpulse enhancements gen- 
erally follow bundles of magnetic field lines different from 
those of enhancements corresponding to the profile com- 
ponents (Fig. |l|- right-hand side). On the contrary, within 
the patchy model both subpulse and profile component en- 
hancements follow approximately the same bundles of field 
lines (Fig. [l] - left-hand side). This should be clear from 
the composite Fig. [l| In fact, subpulse enhancements are 
associated in both models with small subpulse spots S 
following a narrow bundle of dipolar field lines, while the 
profile components correspond to the conical structures 
C in the conal model, and to the narrow patches PA en- 
closing subpulse spots in the patchy model. Because of 
the diverging nature of dipolar field lines controlling the 
plasma flow, all discussed emission features S (spots) , SP 
(sub-patches), PA (patches) and C (cones) are frequency 
dependent, and two frequencies v\ and v-i are marked for 
each feature in Fig. |l|. 
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Fig. 2. Simulation of the distribution of separations Aip s 
between subpulse peaks as compared with separation Aip p 
between the component peaks (vertical solid lines) mea- 
sured at two frequencies (400 and 1400 MHz) within the 
framework of the conal model (panel a) and various ver- 
sions of patchy model (panels b - f). The number of pulses 
(out of the total number of 2000 simulated) is shown on 
the vertical axis and the separation in degrees of longitude 
is shown on the horizontal axis. The dashed vertical line 
going through panels b - f is shown to refer to the sepa- 
ration A(p p — 2°. 875 in the pure conal model (panel a). 



2.1. Conal model of pulsar beams 

The frequency dependence of pulsar emission patterns 
within the angular beaming model of subpulse emission, 
and the related conal model for a mean pulsar beam 
is illustrated in the right-hand side of the composite 
Fig. |l|. The subpulse-associated sub-beams corresponding 
to HPBW (half power beam width) of subpulse emission 
(thin small circles), which are called spots and marked 
by S(v), perform a more or less organized circumferen- 
tial motion around the magnetic axis m, in which the 
azimuthal angle o varies with time, while the opening an- 
gle p[y) remains constant. On average, this motion de- 
termines the cones (thick large circles) of the maximum 
mean intensity with the opening angles p\ and pi at fre- 
quencies v\ and V2, respectively. The frequency-dependent 
longitude <fp(y) of the profile component is determined by 
the intersection of the line-of-sight with the average cone 
at the frequency-dependent opening angle p{v). On the 
other hand, the frequency-dependent longitude (p s {y) of a 
subpulse peak is determined by the local maximum inten- 
sity along the cut of the line-of-sight through the subpulse 
spot S(v). When the frequency changes from v\ to V2, the 
subpulse peak longitude changes from ip a {v\) to ip a (v2), 
while the corresponding profile peak longitude changes 
from <fp{v\) to (p p {v2)- Note that the frequency separation 
Aip p = ippiy-2) — ip p (vi) of profile peaks is generally larger 
than the frequency separation Aip s = <p B {y2) — fsi^i) of 
subpulse peaks. This is consistent with observations pub- 
lished by Izvekova et al. (1993| ), confirmed recently by 



GGGK. As we argue in the next section, within the patchy 
beam model Atp p « Aip s , which is quite different from the 
conal model in which Aip p > Atp s . 

We now attempt to quantify our qualitative statements 
presented above. We use the geometrical method of trans- 
ferring the instantaneous intensity distribution from the 
polar cap to the radio emission region along the dipo- 
lar field lines. To be able to compare the results of sim- 
ulations with the observational data of PSR B0329+54 
(GGGK) we adopt the following parameters P = 0.714 s, 
P = 2-10~ 15 , a = 20°, f3 = 4°, v x =0.4 GHz and v 2 = 1.4 
GHz. The conal model is represented by an arrangement 
of 12 equisized and equidistant sparks, each having the 
HPBW diameter T> ~ 0.2r p , circulating around the mag- 
netic pole at a distance of about 2 /3 of the polar cap radius 
r p . The results of simulations are presented in Fig. 2 (panel 
a). The vertical solid line at about 2.9° represents the sep- 
aration A(p p between the peaks of the average component 
measured at the two frequencies. The distribution of the 
separations Aip s of subpulse peaks measured at the two 
frequencies has a width of about 0.8° and peaks around 
1.1°. The apparent difference between the peak of sub- 
pulse separation and the component peak separation i.e. 
Aifp — At/?P cak ~ 1.8°, as well as the width and somewhat 
skewed shape of the distribution, resemble the observa- 
tional data quite well (see Fig. 4 in GGGK, keeping in mind 
that zero position bin in their histograms corresponds to 
cases when Aip p = Aip s ). We have taken into account only 
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subpulses that corresponded to the same spark at both 
frequencies. In terms of the observational data this corre- 
sponds to subpulses correlating at two frequencies (refer 
to the CCF technique used in GGGK). 

2.2. Patchy model of pulsar beams 

Within the patchy model of pulsar beams (LM88; 
Manchester, 1995] ) the observed pulse profile is the prod- 
uct of a "source function" and a "window function" , 
which can be related to the energy/density distribution of 
plasma beams along different bundles of dipolar field lines 
and to properties of the emission mechanism, respectively. 
This model is illustrated schematically on the left-hand 
side of the composite Fig. 0. The subpulse-associated spots 
marked by SP{v) corresponding to the HPBW of subpulse 
emission (thin small circles), occur within the limited ar- 
eas PA(u) called patches (the simplest and probably un- 
realistic model is when a patch PA can accommodate 
just one sub-patch SP). The occurrence of SP(v) within 
PA(y) can be completely random, or more or less orga- 
nized (including drifting). The dependence on frequency 

v reflects the diverging nature of dipolar field lines, thus 

vi < v\ . The bundle of field lines associated with a patch 
PA(y) is only slightly larger than each bundle associated 
with SP(y). This means that enhancements corresponding 
to subpulses and profile components follow approximately 
the same bundles of field lines. Thus, the frequency depen- 
dence of emission patterns in the patchy model (left-hand 
side of Fig. |l|) should be different from that of the conal 
model (right-hand side of Fig. |l|), in which subpulse en- 
hancements (spots) and profile components (cones) gen- 
erally follow quite different bundles of field lines. If the 
patch is relatively small as compared to the entire po- 
lar cap, then the frequency separation of the profile com- 
ponents Atp p — <p p (v2) — fpiyi) an d that of subpulses 
A(p s — tp s {v2) — (Ps{vi) should be about the same. We 
demonstrate this below by means of geometrical simula- 
tions. 

We have calculated a sequence of single pulses and 
average emission again for the case of PSR B0329+54, 
assuming that the subpatch SP(v) is comparable in size 
with the spark associated emission in the conal model, and 
that the patch PA(v) is twice larger than SP{v). Thus, 
the projection of a patch onto the polar cap has a char- 
acteristic dimension P ~ TD ~ 0Ar p . Such a patch en- 
compasses about 30° in magnetic azimuth (roughly corre- 
sponding to the scale presented in Fig.l (left-hand side)). 
The results of simulations are presented in Fig. 2 (panel b). 
The solid vertical line at about 1°.25 represents the sep- 
aration of the profile component associated with PA(v) 
measured at the two frequencies. The narrow distribution 
of subpulse peak separations peaks exactly at the same 
value (we have again taken into account only subpulses 
associated with the same subpatches SP{v) at both fre- 
quencies). It is worth noting that if we chose the sizes of 
SP(v) and PA{y) to be equal (which is the simplest and 



an unrealistic model of the patchy emission) , the distribu- 
tion would be represented by a delta function coinciding 
with the solid vertical line. The case presented in panel 
(b) is inconsistent with the frequency dependence of pul- 
sar radiation patterns (GGGK and reference therein) . 

We now start to increase the azimuthal dimension of 
the patch PA(v) along the cone, keeping the radial dimen- 
sion the same (about 2 spark diameters V or about 0Ar p ). 
Panels c - f in Fig. 2 correspond to the elongation factor 
1.3, 1.5, 2.0 and 2.5, respectively. Thus, the maximum 
elongation is equivalent to the length scale comparable 
with the polar cap radius r p . As one can easily notice from 
panels c - f , increasing elongation results in two effects: (1) 
the separation of the component peaks increases towards 
the value corresponding to the conal model (dashed ver- 
tical line) i.e A(p p — > 3°, while the distribution of sub- 
pulse peak separations Aip s gets broader and broader and 
peaks at correspondingly larger and larger distances from 
Aip p . Moreover, the skewed shape of the distribution be- 
comes more and more apparent in panels (d) and (e), to 
such an extent that the case presented in panel (f) seems 
almost undistinguishable from the pure conal case pre- 
sented in panel (a), except that the width of the distri- 
bution is too large compared with observations (see Fig. 
4 in GGGK). Further increasing of the elongation factor 
beyond 2.5 (corresponding to 60° of magnetic azimuth or 
1 /6 of the full cone) does not practically change the results 
of the simulations. This means that the unrealistic patches 
elongated to a large extent along circles centered on the 
magnetic axis would resemble the conal model. Although 
this conclusion seems trivial, it allows us to constrain some 
characteristics of both patchy and conal models. This is 
described in the two following paragraphs. 

First we can ask: what is the probability that an ade- 
quately large and favourable patch elongated along a cone 
and corresponding to panel (f) mimics the conal model 
presented in panel (a) of Fig. 2. Of course, we have to take 
into account that another similar patch is required on the 
opposite side of the fiducial plane (see Fig.l) to account 
for the second outermost component of PSR B0329+54. 
Let us assume for simplicity that our elongated patch is 
a rectangular figure with the shorter side A 0Ar p and 
the longer side B ~ 1.0r p (see above for an estimate of 
the dimensions). Thus, the surface area of such a patch 
S P atch ~ A-B = 0Ar p and the probability of its occurence 
in any location of the polar cap is V\ = S patc h/ S cap = 
OAr^/ '(irr p ) 2 w 0.13. The probability of occurence of an- 
other such patch somewhere on the remaining part of the 
polar cap is V 2 = S pa tch/(S C ap - S pa tch) ~ 0.15. To es- 
timate the probability of the proper alignment along the 
cone, we can calculate a number of different independent 
orientations of our rectangular figure inscribed into a cir- 
cle of the diameter approximately equal to A = r p . One 
can easily show that the number of independent orienta- 
tions is approximaly 27r(r p /2)/(0.4r p ) ~ 8 and thus the 
probability of alignment along a cone Vo ~ 0.12. So far, 
the resultant probability is V = V x ■ T 2 ■ Vq ~ 3 ' 10 ~ 4 - 
The requirement of having two such patches symmetricaly 



6 



J. Kijak & J. Gil: Structure of pulsar beams: conal versus patchy 



placed with respect to the fiducial plane only decreases 
the final probability. Thus, we can conclude that mimick- 
ing the conal radiation pattern by a specially arranged 
patchy distribution is extremly unlikely, at least in the 
case of PSR B0329+54 and other pulsars showing similar 
frequency dependence of emission patterns (GGGK and 
references therein). 

We can now approach the results of our elongation 
exercise from a different angle. Since further elongation 
beyond the case presented in panel (f ) does not change the 
obtained distribution, we can constrain the possible spread 
of the cone in the radial direction as compared with the 



malized polar cap for short and long period pulsars, 
respectively. 

(ii) A similar problem concerns the putative conal struc- 
ture of pulsar beams which HM01 did not reject a pri- 
ori. However, they assumed that adding the normal- 
ized polar caps of different pulsars would not smear the 
conal structure, if it existed. It seems that the conal 
patterns are almost certainly different in different pul- 
sars. MD99 showed that the pulsar emission beams fol- 
low a nested cone structure with up to three distinct 
cones, although only one or more of the cones may be 
active in a given pulsar (see also Rankin, 1993 ; Gil ct 



ideal ( and probably not realistic) case presented in panel 
(a). Remembering that the initial patch size PA (panel 
b) was twice SP, with SP corresponding to the spark 
size V, we conclude that the realistic version of the cone 
can accommodate up to two sparks in radial dimension. 
This means that a locus of maximum intensity within the 
average pulsar beam has the form of a narrow ring rather 
than a circle. 

2.3. Beam reconstruction techniques: 

Recently Han & Manchester (2001; HM01 hereafter) at- 
tempted to reveal the shape of pulsar radio beams. They 
claim to have constructed a two-dimensional image of the 
"average" mean pulsar beam using a special technique ap- 
plied to all available multi-component pulse profiles with 
good quality polarimetric data. They mapped the ob- 
served profile intensity onto the line-of-sight crossing the 
normalized polar cap at the normalized impact angle (3 n 
estimated from the polarization angle swing. To include a 
second dimension (perpendicular to the line-of-sight), they 
broadened the distribution in latitude applying a Gaussian 
to each longitudinal sample. Adding all pulsars together 
they obtained a global average beam pattern, which (after 
normalization to correct for the nonuniform distribution of 
P n ) they believe represents the global mean pulsar beam. 
HM01 concluded that their results are consistent with the 
patchy rather than conal beam model (see their Fig. 4). 
However, as we demonstrate in items (i)-(iv) below, their 
beam reconstruction technique is not general enough to 
reveal the true structure of pulsar beams. 



(i) HM01 projected all emission features onto a normal- 
ized polar cap, ignoring the dependence of the beaming 
angles on the emission altitude, which most probably 
depends both on the radio frequency and on the pul- 
sar period and its derivative (Kijak & Gil, 1997, 1998; 
Kijak, 2001). Taking into account the diverging na- 



al., 1993] ) . Also GS00 argued that the number of cones 



ture of dipolar field lines with increasing altitudes, the 
projection of the emission pattern onto the polar cap 
must be performed more carefully. HM01 used data 
at frequencies between 600 MHz and 1.6 GHz, which 
probably takes care of frequency dependence of the 



op ening angles p(v) 



,1/2/ 



However, the 



per od dependence p(p) oc p~ u i is stronger and can 



result in a broader spread of projections onto the nor- 



is a function of the basic pulsar parameters P and P. 

(iii) The analysis of HMO 1 is based on the orthogonal nor- 
malized impact parameter (3 n = 890/ P90 , where 890 
and pgo are the impact angle 8 and the beam radius 
p computed for the inclination angle a = 90° (LM88) . 
The orthogonal normalized impact parameter 8 n dif- 
fers from the actual normalized impact parameter 8/ p 
in many aspects (e.g. Gil et al., 1993). We would like to 
emphasize here that [3/ p depends on the pulsar period 
P, the inclination angle a and the observing frequency 
v. These dependences, which can affect the emission 
patterns, are not accounted for in the HM01 analysis 
based on the orthogonal impact parameter. 

(iv) HM01 included the core components in their analy- 
sis. This is another possible source of confusion, since 
the core components are rather randomly placed with 
respect to the midpoint of the overall profile. For this 
reason MD99 excluded the core components from their 
analysis, which revealed a nested cone structure of pul- 
sar beams. 

Given the problems listed above, we conclude that the 
results of the HMOl's beam reconstruction are illusive. 
The lack of an apparent conal structure in their 'global' 
pulsar beam does not exclude the existence of such struc- 
tures in particular pulsars, expecially if they are deter- 
mined by physical and geometrical factors, which may 
vary a great deal among different pulsars. HM01 demon- 
strated that the conal emission is not confined to a single 
region at the beam boundary. They also concluded that 
if multiple cones exist, they are at different radii relative 
to the beam radius in different pulsars. This is exactly 
what the multi- nested cone model of GS00 predicts. Both 
the number of cones and the relative radius of a given 
cone depend on a pulsar period and its derivative in their 
model. 

As an example of pulsar modelling within the multi- 
nested cone scenario (GS00), we use the case of PSR 
J1834— 0426. The patchy model for this pulsar was pre- 
sented by HM01 (their Fig. 1). This pulsar has a very 
broad profile with the pulse width W ~ 130 degrees of lon- 
gitude, which implies a small inclination angle a <C 90°. 
Assuming a real istic emission altitude r$ ~ 40 ( Kijak & 
Gil, 1997| , |1998|) we estimated the observational angles 
a ~ 10° and 8 ~ 2°. 5 (consistent with a = 7.9° and 




Fig. 3. The proposed geometry of observation and the instantenous arrangement of sparks on the polar cap for conal 
model of pulsar beam of PSR J1834— 0426. Notice that about 30% of the line-of-sight (1-of-s) trajectory covers the 
beam, resulting in very broad (~ 130°) profile of this pulsar. 



(3 — 1.6° given by LM88). Figure 3 shows an instan- 
tenous arrangement of sparks on the polar cap of PSR 
J 1834— 0426, obtained from the complexity parameter of 
GS00. The circumferential motion of these sparks results 
in the average structure of two nested cones (e.g. Fig. 6 in 
GGGK). Notice that about 35% of the line-of-sight trajec- 
tory stays within the beam, which leads to a broad ~ 130° 
average profile presented in Fig. 4. 



3. Discussion and conclusions 

In this paper we explore a geometrical method of pul- 
sar radiation simulation, based on two well-justified as- 
sumptions: (i) the elementary coherent radio emission is 
narrow-band, and the emission altitude depends on both 
the frequency and the pulsar period, (ii) the emission is 
relativistically beamed tangently to dipolar magnetic field 
lines. We have considered two competitive models of the 
organization of pulsar emission beams: the conal model, 
in which enhancements related to subpulsc emission in 
single pulses are distributed along the cones correspond- 
ing to maximum average intensity, and the patchy beam 
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Fig. 4. The observed (solid line) and simulated (dashed line) mean profile of PSR J1834— 0426. The conal model 
calculations corresponding to 610 MHz were performed following the instantaneous energy distribution (presented in 
Fig. ||). The geometry of observations was determined by the inclination angle a = 10° and the impact angle /3 — 2°. 5. 
The adopted emission altitude r(610 MHz) « 30R. The 610 MHz observational data ( Gould fc Lyne, 199§| ) were taken 
from the European Pulsar Network data base (http://www.mpifr-bonn.mpg.de/div/pulsar/data/) 



model in which subpulse enhancements corresponding to 
the component of the mean profile are confined to the 
patchy area limited both in azimuthal and radial dimen- 
sions. We examined the consistency of these rival models 
with the variety of observational data. 

We have argued that a number of observational prop- 
erties of pulsar radio emission, namely: (i) binomial dis - 
trib ution of the opening a ngles ( Rankin, 1993| ; Gil et al 



1993; Kramer et al., 1994 ); (ii) high impact angles corre- 



sponding to single and double profile pulsars (MD), and 
(iii) different frequency dependence of a subpulse and cor 
responding profile component longitudes ( Izvekova et al. 
199 3j; pi fe Krawczyk, 1996fc pil et al., 2002j ), strong!^ 



suppor t the conal model ot pulsar beams. The alterna- 
tive patchy beam model is inconsistent with these obser- 
vational properties of pulsar radiation. 

We have also demonstrated that the beam reconstruc- 



is not capable of revealing the true structure of pulsar 
beams. In fact, their formalism assumes implicitely that 
neither the radio emission altitude nor the number of pu- 
tative nested-cones and their locations within the pulsar 
depends on the pulsar period. The lack of an apparent 
conal structure in their "global beam" does not exclude 
the conal beam model. Thus, the results of the HM01 anal- 
ysis provide no strong evidence of patchy beam structure 
in pulsars. 

We tend to favour the version of the conal model in 
which the relationship between the subpulse-associated 
beams and cones of the average emission is established 
through the phenom e non of the E x B drift ([Ruderman 



fc Sutherland, 1975; Dcshpandc fc Rankin, 199E, 2001 



tion te chnique developed by Han fc Manchester (2001 ) 



GS00), which forces the spark filaments of plasma to ro- 
tate slowly around the magnetic axis. This "spark model" 
of radio p ulsars was recently tested statistically by Fan et 
al. (200l[) . They showed, by means of Monte Carlo simu- 
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lations, that various pulsar parameters can be reproduced 
if both the spark dimension and their mutual separation 
are approximately equal to the height h of the polar gap 
(Ruderman & Sutherland, 1975, GSOO), or consequently, 
the maximum number of sparks along the diameter of the 
polar cap with the radius r p is N max ~ r p /h. Thus, their 
conclusions are consistent with the assumptions used in 
this paper. 

We note that evidence of a relationship between 
drifting subpulses and the conal structure of mean pul- 
sar beams already exists in the literature. Hankins & 



Wolszc zan (1987 ) examined three pulsars with triple aver- 
age profiles showing subpulses drifting across the full pulse 
window (including the central component). They found 
that these pulsars are consistent with two nested cones 
of emission, eac h associated with a prominent subpulsc 
drift. Moreover, Proszyhski fc Wolszczan (1986 ) showed 
that in complex profile pulsars there is a strong correla- 
tion between drifting subpulses associated with different 
profile components. Such correlations are natural within 
the E x B induced conal model, but inconsistent with the 
patchy model of pulsar beams. 

Finally, we suggest that more pulsars should be ob- 
served in the single pulse, simultanenous dual frequency 
mode. Our simulation method illustrated and described in 
Fig. 2 can be adapted to the analysis of such data in or- 
der to ultimately discriminate between conal and patchy 
beam models in pulsars. 
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